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Abstract 

The new multi-color BVRI photometric light curves of the short-period eclipsing 
binary GSC 3576-0170 were obtained on two consecutive nights (October 5 and 6, 
2009). With the 2003 version of Wilson-Devinney program, the precise photometric 
solutions are derived for the first time. The result shows that GSC 3576-0170 is a 
semi-detached binary system with a large temperature difference of approximately 
1490 K. The light-curve distortions are further explained by a hot spot on the 
secondary component through mass transfer via a stream hitting the facing surface 
of the secondary component. By analyzing all available light minimum times, we 
also derived an update ephemeris and found for the first time a possible periodic 
oscillation with an amplitude of 0.0038 days and a period of 4.3 years. The periodic 
oscillation could be explained either by the light-time effect due to a presumed third 
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component or by magnetic activity cycle of the system. 

Key words: stars: binaries: close - binaries: eclipsing - stars: individual (GSC 
3576-0170) 97.80.Fk 97.80.Hn 97.20.Jg 



1 Introduction 

GSC 3576-0170(Por6=0'^.405, Gl V) is a new near-contact solar-type eclipsing 
binary, which displays asymmetric light curves (two different light maxima) 
(Nelson et al. 2006). Therefore, it is a very intriguing target for understanding 
the property of the system. 

GSC 3576-0170 was discovered to be variable by Nelson et al. (2006). They 
attempted a preliminary light-curve synthesis with Wilson-Devinney modeling 
and preferred a detached model with a mass ratio in the range 0.15 — 0.35 
and an orbital inclination of 65-70°. For the period variation, they invoked a 
quadratic function to fit the values of the observational times of light minimum 
- calculational times of light minimum {O — C). Subsequently, new minimum 
times of GSC 3576-0170 were published by several astronomers (Nelson et al. 
2006; Hiibscher & Walter 2007; Hiibscher et al. 2006; Brat 2007; etc.). 
In this paper, we present our new V, R, and / LCs of GSC 3576-0170 which 
were analyzed using the 2003 version of Wilson-Devinney code (Wilson & 
Devinney 1971; Wilson 1979, 1990, 1994; Wilson & Van Hamme 2004). We 
accumulated all available times of light minimum and discussed the period 
change. 

Email address: liyjzhangOhotmail . com, phone: +86 851 362 7662, fax: 
+86 851 362 7662 (Liyun Zhang^'^). 
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Table 1 



Minimum times of GSC 3576-0170. 
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References - (1) Nelson et al. 2006; (2) Hiibscher & Walter 2007; (3) Brat 2007; (4) 
Hiibscher et al. 2006; (5) Present paper. 
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2 New CCD photometric observations of GSC 3576-0170 



New CCD observations were made on October 4 and 5, 2009 with the 85 cm 
telescope at Xinglong Station of National Astronomical Observatories of China 
(NAOC). The telescope was equipped with a standard Johnson-Cousin-Bessel 
multi-color CCD photometric system built on the primary focus (Zhou et al. 
2009). The camera was equipped with a 1024 x 1024 CCD and the field of 
view was about 16'. 5 x 16'. 5 (Yang et al. 2009). All CCD images were reduced 
by means of IRAF package in the standard fashion. The star (q;2ooo = 20 : 
23 : 44.56; ^2000 = 46 : 51 : 33.9) and GSC 3576-0970 (ctaooo = 20 : 24 : 04.04; 
^2000 = 46 : 53 : 21.8) were chosen near the target as comparison and check 
stars, respectively. The magnitude of these stars were determined using the 
Daophot subpackage of IRAF and the errors were found to be lower than 0.01 
mag in all bands. Corresponding light curves are displayed in figure 1, where 
the phases of data points are calculated using the new light ephemeris (the 
equation 1). 



3 Orbital Period Variation for GSC 3576-0170 

From our observations, a new secondary minimum (2455109. 0636±0. 0025) and 
a primary minimum (24551 10.0739± 0.0015) for GSC 3576-0170 are obtained 
by means of parabola fitting. To discuss periodic variation, we compiled all 
available light minima, which are listed in table 1. Using the linear least- 
squares method, an update linear ephemeris is given by equation(l): 

Min.I = JD2452795.8754(±0.0007) + 0''.4050051(±0.0000004)£; (1) 
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Fig. 1. CCD photometric observations of CSC 3576-0170. 
Based on our new ephemeris, the — C values of minimum times are also listed 
in table 1. The calculation and the plot are provided in figure 2. As can be 
seen from figure 2, a cyclic variation might be present. To clarify the period 
change, we used a sinusoidal function for fitting which led to the following 
equation: 

0-C = 0'^.0038(±0.0006) x sin(0°.093(±0°.005) x E - 27°.1(±1°.3)X.2) 

The fit is also plotted in figure 2 with solid line. The sinusoidal term of equa- 
tion (2) reveals a periodic oscillation with an amplitude oi A = C^.OOSS. Using 
T=360°xP/(X', where P denotes the orbital period of GSC 3576-0170 in years 
and u is the coefficient of E; the period of the oscillation T is determined to 
be 4.3 ±0.2 years. The size of the periodic variation AP/P=1.5xlO^^ can be 
computed in days using the equation AP/ P = Ax \J2[\ — cos{2ti x P/T)]/P 
(Rovithis-Livaniou et al. 2000), where A is the semi-amplitude of the light- 
time effect, and T is the period of the light-time effect. 
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4 Photometric analysis 



For our light curves of GSC 3576-0170, two light maxima are almost equal in 
BVRI bands (the differences are all less than 0.01 mag). This implies that it 
can be used to determine reliable photometric parameters. During the solution 
process, the four light curves are simultaneously analyzed. The preliminary 
values of these orbital parameters are taken from the preliminary photometric 
solutions derived by Nelson et al. (2006). The details of the procedure of pho- 
tometric solution are similar to those of our previous work of RT And (Zhang 
& Gu 2007), DV Psc (Zhang et al. 2010), and KQ Gem (zhang 2010). 
To be on the safer side, we explored a wide range based on the estimation 
of Nelson et al. (2006). The mass ratios were fixed at 0.2, 0.25, 0.3, 0.35, 
0.4, 0.45, 0.5, 0.6, and 0.7, respectively. For each assumed value of q. the 
light curves were solved based on mode 2 (detached mode) of Wilson and 
Devinney's DC Program. According to the spectral type Gl V of the pri- 
mary (Nelson et al. 2006) and the spectra-effective temperature relation of 
Cox (2000), the effective temperature of the primary Ti was fixed at 5865. 
During the solution process, we assumed synchronous rotation and zero ec- 
centricity. Simple treatment was used to compute the reflect effect, and the 
linear limb-darkening law was employed to compute the limb-darkening ef- 
fect. The bolometric albedo ^41=742=0. 5 (Rucindki 1973), the hmb-darkening 
coefficients xib = 0.742 X2b = 0.906, xiy = 0.608 X2v = 0.763,a;iij = 0.503 
X2R — 0.629, Xii — 0.413 X21 — 0.512 (Van Hamme 1993), and the gravity- 
darkening coefficients gi — 92 — 0.32 (Lucy 1967) are set for the primary and 
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the secondary in the usual manner. The solutions from several runs were ob- 
tained for all assumed mass ratios. The relation of the sum of weighted square 
deviation (J2i{0 — C)^) with mass ratio q is illustrated in figure 3, where the 
lowest value of was found at q = 0.45. Therefore, it indicates that the most 
likely mass ratio appears to be approximately 0.45. The photometric solutions 
for 0.45 are listed in table 2 and corresponding light curves are plotted with 
dashed fines in figure 4. 

Comparing the observed LCs with the theoretical ones for the likely mass 
ratio (see figure 4) shows enhancement (brighter) on the shoulders of the sec- 
ondary eclipse. This kind of enhancement may be caused by mass transfer 
from the primary component. To account for this phenomenon, we adopted 
a semi-detached model with the addition of a hot spot on the secondary. For 
this purpose, we changed the mode from 2 to 4 (semi-detached mode with 
lobe-filling primary), and continued the DC calculation until the final pho- 
tometric solution of was derived. At the same time, the mass ratio q were 
released as a free parameter and continue the DC calculation based on our 
initial solutions of mode 2. To avoid the correlations among the adjusted pa- 
rameters, they were divided into two subsets: the orbital parameters and spot 
parameters. For the spot parameters, the spot latitude is assumed to be 90°, 
which means that spot center is on the equator of the component. The spot 
temperature is assumed to be the same as that of the primary (5860 K) due 
to the mass transfer, that is transferred from the primary component, would 
be about that of the primary component. Therefore, we only adjusted the re- 
maining two parameters of the spot, the spot longitude longitude spot ^-nd the 
radius radiusgpot- The solutions from several runs were derived. Correspond- 
ing photometric solutions are listed in table 2. The theoretical light curves are 
plotted in figure 4 with solid lines, and corresponding configurations of CSC 
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Table 2 

The results of LC analysis for GSC 3576-0170. 



Element 


detached 


semi-detached 




5865K a 


5865K a 


q 


0.450a 


0.467±0.003 


i 


64. 32ili0.69 


64. 63±0.14 




4369ib46 K 


4375±27 




2.702±0.004 


2.798a 




2.880±0.023 


2.790±0.005 


Lib/{Li + L2)b 


0.9587±0.0002 


0.9444±0.0002 


Livl {Li + L2)v 


0.9302±0.0004 


0.9073±0.0004 


Lir/{Li + L2)r 


0.9088±0.0005 


0.8807±0.0005 


Lii/{L^ + L2)i 


0.8830=b0.0007 


0.8485±0.0006 


ri(pole) 


0.4368±0.0008 


0.4200±0.0005 


ri(side) 


0.4684±0.0010 


0.4466±0.0005 


ri (back) 


0.5015±0.0014 


r\ AfJ A ^ 1 ( \ f\C\f\ r" 

0.4743±0.0005 


r2(pole) 


0.2743±0.0037 


0.2984±0.0013 


r2 (side) 


0.2840±0.0043 


0.3119±0.0016 


r2 (back) 


0.3073±0.0060 


0.3470±0.0027 


Latitudegpot 




90° a 


Longitudcspot 




181.°9±0.°2 


RadiuSgpot 




27.°8±0.°5 


Temperature spot 




5865 K a 


no - c)i 


0.7188 


0.2694 



Parameters not adjusted in the solution are denoted by a mark "a" 



3576-0170 in phases 0.0, 0.25, 0.5 and 0.75 are shown in figure 
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Fig. 2. The 0-C diagram for the minimun times of GSC 3576-0170. The solid hne 
represents the sin fitting. 
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Fig. 3. X] - q curves of GSC 3576-0170 
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1.4 



theoretical light curves for detached mode 
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Fig. 4. The observational and theoretical light curves of GSC 3576-0170. The circles 
represent the observational data. The dashed and solid lines represent the theoretical 
light curves for detached and semi-detached mode, respectively. 
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Fig. 5. The configurations of GSC 3576-0170 in phases 0.0, 0.25, 0.5 and 0.75. 
5 Discussions and Conclusions 



Our new LCs and the period variation of GSC 3576-0170 are included in 
present investigation. 
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5.1 Cyclic magnetic activity or the light-time effect due to the third Body? 

For GSC 3576-0170, there seems to be the sinusoidal oscillation of the O- 
C residuals. The oscillating characteristic may be caused by the hght-time 
effect due to the existence of the third body orbiting the eclipsing binary or a 
possible result of magnetic activity cycles of the system. 
On one hand, the period oscillation may be explained by the light-time effect. 
Assuming that the orbital of the presumed third body is circular, we can 
obtain the mass function for the third body f(M3)=0.016(±0.001) M© using 
the following equation: 

./X 47?^ . (Masinio)^ , , 

where Mi, M2 and M3 are the masses of the binary, and the third body, re- 
spectively. The mass of the third component can also be calculated with the 
above equation, which depends on the orbital inclination. The minimal mass 
M^^rnvn IS 0.38 M© wheu i'=90°. 

On the other hand, the period oscillation may be accounted by magnetic cycle 
of the system (Applegate 1992; Lanza et al. 1998; Lanza & Rodono 1999). The 
orbital period change corresponding to a variation of the quadrupole moment 
is given by AQ = -(Ap/p)x(Ma^/9) (Applegate 1992) where a is the semi- 
major axis of the binary orbit and M is the mass of the active component. 
Since no spectroscopic solutions are available in literature GSC 3576-0170, its 
absolute parameters can not be directly determined. We estimated the pri- 
mary mass 1.024 as M©, and the radius as 1.06 R© by assuming the primary 
component to be a normal and main-sequence GIV star. Based on our pho- 
tometric solutions, the mass of the secondary component is M2 — 0.478 M© 
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and the separation between the two components is a = 2.37 R©. Therefore, 
the quadrupole moment AQ is calculated to be 0.5x10^° g cm^ for the pri- 
mary. This value is similar to the typical values for active binaries (Lanza & 
Rodono 1999). Therefore, the magnetic activity cycle is a possible mechanism 
to explain the variation of period. 

5.2 Long-term change and their evolutionary status 

For our observations, the two light maxima are basically equal. However, the 
light curves show that the first quadrature is brighter than the second one in 
2003 (Nelson et al. 2006). These indicate a variability on a time scale of about 
six years. This variation might be attributed to surface magnetic activity of 
one or both stars. More importantly, GSC 3576-0170 is perhaps associated 
with possibly thick deep convection according to their spectral types Gl V. 
Thus, it confirms that magnetic activities may produce the sinusoidal varia- 
tion of the period. 

Our result showed that GSC 3576-0170 is a marginal-contact binary system 
with a large temperature difference of about 1490 K between the two compo- 
nents, which belongs to the subclass of VIOIO Oph binaries (Shaw 1994). It is 
similar to several semi-detached systems with a possible third body, such as 
GW Tau (Zhu & Qian 2006). They may be potential candidates of binaries 
that will become overcontacted or remain in the broken- contact in accordance 
with the thermal relaxation oscillation theory of overcontact binaries (Lucy 
1976; Flannery 1976; Robertson & Eggleton 1977; etc.). 
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Of course, our solutions are based on photometric observations only. For better 
understanding of the properties and the evolutionary state of GSC 3576-0170, 
spectroscopic observations are needed. Moreover, the 0-C residual of GSC 
3576-0170 is about seven years and there is an observational gap. Therefore, 
it is too early to decide the character of the periodic variation, and it might 
take another 10-20 years are needed for confirmation. 
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